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Abstract 

This paper considers the problem of estimating the delays of a weighted superposition of pulses, called 
stream of pulses, in a noisy environment. We show that the delays can be estimated using a tractable 
convex optimization problem with a localization error proportional to the square root of the noise level. 
Furthermore, all false detections produced by the algorithm have small amplitudes. Numerical and in- 
vitro ultrasound experiments corroborate the theoretical results and demonstrate their applicability for 
the ultrasound imaging signal processing. 


1 Introduction 

In many engineering and scientific problems, we acquire data that can be modeled as a weighted super¬ 
position of pulses (kernels), and aim to decompose it into its building blocks, frequently called atoms. Typical 
examples are ultrasound imaging |331 , 531 and radar [3], where the measurements are echoes of the emitted 
pulse, that are reflected from different targets. When attempting to decompose the stream of pulses there 
are two main concerns: the robustness of the estimation and its degree of localization. While the localization 
of the estimation determines how close is the estimation to the real location of the targets, the robustness of 
the estimation determines the degree to which reliable reconstruction can be performed under different noise 
conditions. This paper investigates the localization of the decomposition of stream of pulses with application 
to ultrasound imaging. 

Ultrasound imaging is currently one of the most widely-used medical imaging modalities worldwide. 
The popularity of this modality can be attributed to a variety of appealing characteristics, such as cost- 
effectiveness, portability and being practically harmless. In the past few years there is a shift in ultrasound 
systems from hardware based beam-forming and signal processing to software based implementations. In 
these new systems the RF data is accessible with greater ease, enabling implementation of advanced signal 
processing algorithms. Although ultrasound scans are sufficiently accurate for many applications, the resolu¬ 
tion of this modality is limited by the finite bandwidth of ultrasound probes and the aperture of the transducer 
which determines the axial and lateral resolution, respectively. The limited resolution of ultrasound scans 
obscures important details and limits the clinical useability of the modality. Another inherent feature of ul¬ 
trasound imaging, being a coherent imaging modality, is speckle noise. Biological tissues are characterized by 
numerous structures much smaller than the wavelength of the ultrasonic pulse. These structures induce local 

1 These two authors contributed equally to this work. 


1 



inhomogeneity of the acoustic impedance. As a result, the ultrasonic pulse is reflected from many indepen¬ 
dent scatterers in each resolution cell, defined by the -6dB of the main beam in the axial, lateral and out-of 
plane dimensions. The combination of these reflections creates a complicated constructive and destructive 
interference pattern, named speckle noise. Even though the small reflectors creating the speckle noise cannot 
be resolved, the localization of strong reflectors is important for many clinical application, motivating the 
development of deconvolution and super-localization methods [261 ;T[ [25]. 

The ability of different biological structures to reflect the ultrasonic wave is commonly represented by a 
3D function, called the reflectivity function. When assuming weak scattering and linear propagation, the 
received signal (RF signal) can be modeled as a convolution of the reflectivity function and the point spread 
function (PSF) of the ultrasound scanner [24j . Even in cases where non-linear propagation is dominant, 
the convolution model can be assumed to hold for short signal sections. The automatic extraction of such 
sections from the received RF data was discussed in several papers (e.g. [251124] ). Even though the ultrasonic 
signal results from a convolution of the 3D PSF and the 3D reflectivity function, the deconvolution problem 
is usually solved by reducing the dimension of the problem and estimating the hypothetical ID reflectivity 
function from each scan line, called a-line. 

We consider a signal (scan line) which consists of a stream of pulses, i.e. 

y [k] = ^2 c m g a [k - k m ] + rj[k], k€ Z, c m e R, (1.1) 

771 

where g a [fc] := g (yyr), g(t) is an admissible kernel, a notion which is defined in Section 2 1/N is the uniform 
sampling interval, a > 0 is a scaling parameter and g is additive noise with \\g\\ := ^ • We aim to 

estimate the true support K := {k m } and the weights {c m } from the measurements {y [k ]}. This paper does 
not deal with sampling techniques and any error due to the sampling process (e.g. aliasing) can be treated 
as noise. 

The model <H3 can be represented as a convolution model 


y [k] = {9c 7 * x) [k] + g[k], 

where V represents a discrete convolution, and 

(1.2) 

x [/c] — ^ ^ [/c km ], 

(1.3) 
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where <5 [ k\ denotes a Kronecker delta function. Our results can be applied directly to deconvolve richer 
classes of signals, as elaborated in Section [2] 

A well-known approach to decompose the signal into its atoms is by using parametric methods such as 
Prony, MUSIC, Matrix Pencil and ESPRIT pTHI [25] [22] [25] , However, the stability of these methods is not 
well-understood, although some steps towards the study of the non-asymptotic behavior have been taken 
recently ES] (27]. An alternative way is to exploit standard compressed sensing and sparse representations 
theorems, relying on the sparsity of the signal [min]- These methods require dictionaries with low coherence, 
and hence highly limiting the resolution of the solution. 

Recently, it was shown that under a kernel-dependent separation condition on the support K (see Defi¬ 
nition 2.11 a signal of the form (1.31 can be estimated robustly by solving an l\ minimization problem J7 . 
This fundamental result is presented in Theorem 2.5 This paper builds on the results of |7] and follows 


m to derive new results, revealing that the recovery of the i\ minimization is well-localized. That is to 
say, the support of the recovered signal is clustered around the support of the underlying signal (1.31 as 
presented in Theorem |2.10| The proof relies on the existence of a special interpolating function, constructed 
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as a super-position of the kernel and its derivatives. Similar techniques were developed in the field of sparse 
recovery, and specifically as the pillars of the super-resolution theory nu na Sana is EE lama. 

As part of the increasing interest in the field, several previous works explored the localization properties 
of super-resolution problems. In Ha ia an, the authors focused their attention merely on the special cases 
of low-passed convolution kernels in Fourier and algebraic polynomials spaces. In this work we treat a much 
broader family of convolution kernels, named admissible kernels (see Definition 2.2). The thorough works 
Hang investigated the localization behavior when the noise level drops to zero. In contrast, the analysis 
in this paper, as presented in Theorem |2.10| holds for all noise levels. Before proceeding, we would like to 
mention that support recovery of sparse signals was studied in many different settings and approaches. For 
instance, isaia analyzed the support recovery of compressed sensing problems using an information-theory 
point of view. 

We validate our theoretical results using ultrasound imaging experiments. These experiments include 
both numerical simulations and in-vitro phantom scans. The numerical simulations enable the validation 
of the main result presented in this paper using numerous reflectivity functions with varying spacings and 
amplitudes. On the other hand, the phantom scans were used to test this method under the conditions of real 
ultrasound signals and a reflectivity function with known geometry. Together these experiments corroborate 
the theoretical results and show their relevance to ultrasound signal deconvolution. 

The outline of this paper is as follows. In section [2j we recall the results of [T1 and present our main 
theoretical result. Section [3] is devoted to the proof of the localization result. Section [3] describes the setup of 
the experiments we have conducted and presents the results, and Section [5] concludes the work and outlines 
some possible future research directions. 


2 Main Result 


In a previous work [7], it was shown that robust recovery of a signal (1.31 from a stream of pulses (1.21 is 
possible if the support K is sufficiently separated and the convolution kernel g(t) meets some mild localization 
properties. To this end, we recall two interrelated fundamental definitions. A signal that its support does not 
satisfy Definition 2.1 cannot be decompose stably . For instance, a signal of the form of y [k] = g a [k — k\] — 
g a [k — k^] where k\ and &2 are adjacent locations on the grid will be drown in a miniscule noise level unless 
the kernel g a is very localized. Definition |2.2| formulates the localization properties of the kernel. Roughly 


speaking, it requires that the kernel and its first three derivatives will decay sufficiently fast and that it is 
concave near the origin. 


Definition 2.1. A set of points K C Z is said to satisfy the minimal separation condition for a kernel 
dependent v > 0 and given scaling parameter a > 0 and sampling spacing 1/iV > 0 if 


min I ki — kA > Nva. 

Definition 2.2. A kernel g is admissible if it has the following properties: 

1. g £ C 3 (R), is real and even. 

2. Global property: There exist constants Cn > 0,£ = 0,1,2,3 such that |g^(t)| < C#/ (l + f 2 ), where 
gW (t) denotes the I th derivative of g. 

3. Local property: There exist constants 0 < e, /J < v such that 
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Kernels 

D 2 

Gaussian := e 2 

Cauchy: = 

Co 

1.22 

1 

c\ 

1.59 

1 

c 2 

2.04 

2 

C 3 

2.6 

5.22 

3 (2) ( 0) 

-1 

-2 

empirical v 

1.1 

0.45 


Table 2.1: The table presents the numerical constants of the global property in Definition \2. 2\ for the Gaussian 
and Cauchy kernels. Additionally, we evaluated by numerical experiments the minimal empirical value of u, 
the separation constant of Definition \ 2. l\ for each kernel. 


(a) g(t) > 0 for all \t\ < e and g(t) < g(e) for all |£| > e. 

(b) g^ 2 \t) < —j3 for all |£| < £. 

_(2 

Remark 2.3. Two prime examples for admissible kernels are the Gaussian kernel g(t) = e~ "3~ and Cauchy 
kernel g(t) = The numerical constants associated with those kernel are presented in Table 2.1 

can be somewhat weakened to \gW (f)| < Ct/{ 1 + |£| 1+s ) 


2.2 


Remark 2.4. The global property in Definition 

for some s > 0. In this case, the separation condition would become dependent on s. 

The kernel-dependent separation constant v is defined as the minimal constant guaranteeing the existence 


of some interpolating function as described in Section 3.1 of [7j (see Lemma 3.1). The minimal separation 
constant for the Cauchy and Gaussian kernels was evaluated numerically to be 0.45 and 1.1 (see Table 2.11. 
In this paper, we consider recovery by solving the following tractable convex optimization problem: 


nun 

XG4 (Z' 


PIL subject to || y-(g a *x) 


\,<s. 


Denote the solution of (2.11 as 


i[k] = y ^c n S 


k — k r 


K := 


■= jfcn} • 


( 2 . 1 ) 


( 2 . 2 ) 


We are ready now to state the main result of |7|- This result shows that the recovery error of (2.1) is 
proportional to noise level. 


Theorem 2.5. Let y be of the form of 0 and let g be an admissible kernel as defined in Definition |^.^| 


If the signal’s support K satisfies the separation condition of Definition 2.1 for a, N >0, then the solution x 

167 % 


of (2.1) satisfies (for sufficiently large v) 


Lr — £ < 


P 


where 


7 = max {TVcr, e 1 } . 


Corollary 2.6. In the noiseless case 


6 = 0, (2.1) results in exact recovery. 


(2.3) 
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Remark 2.7. A tighter bound on the recovery error is presented in Theorem 2.12 in |7j. For the sake of 
clarity, we presented here a simplified version of the stability result. 

Remark 2.8. We chose to focus on the univariate results as we aim to apply it on ID ultrasound scans. 
However, we stress that an equivalent result holds for the bivariate case. Namely, the recovery error is 
proportional to the noise level <5 and depends on the kernel localization properties. 

For a fix k m £ K, let 

Z ra := {k £ Z : | k — k m \ < Nea} . (2-4) 

Accordingly, we define the partition of Z to disjoint subsets Z = Z near U Z f ar , where 

^near := ^{m: ( 2 - 5 ) 

^far : = ^near 

= {k £ Z : \k — k m | > Neaykm £ K} . 


As the notation implies and since e < v, the set T, near is composed of the grid points that are close to a 
support location, whereas Z f ar is composed of all grid points that are located sufficiently far from all support 
locations. 

This paper focuses on the localization properties of the solution of (2.11. Therefore, we present the 


following important corollary which states that the superfluous spikes, far away from the sought signal, have 
low amplitudes. 


Corollary 2.9. Let us define the set Kf ar '■= Zfar H K, where T,f a 
respectively. Under the conditions of Theorem\2.5\ we have 


and K are given in (2.5) and (2.2), 




E 

■k n eK fa 


< 


167 2 


( 5 . 


■} 


Proof. By definition, x[k\ =0 for all k £ Kfar- Therefore, using Theorem 2.5 we get 


E 


E 


{n: k n (E.Kfar\ \n-kn€Kfar\ 

< J2\x[k]-x[k]\ 


< 


fcez 

I67 2 


( 5 . 


□ 


Although the optimization problem ( |2.1[ ) is defined on a discrete grid, it was proven that it converges to a 
solution on the continuum (in the sense of measures) as the discretization becomes finer [31!. Additionally, the 
behavior of the discrete optimization problem solution, when the underlying signal is defined on a compact 
domain is analyzed in detail in mm- 

We want to emphasize that our model can be extended to other types of underlying signals, not necessarily 
a spike train as in (1.3). For instance, suppose that the underlying signal itself is a stream of pulses, 


namely a signal of the form x[k] = c m ga 2 [fc ~ k m \. In this case, the measurements are given as y [£;] = 
(g a 1 * gcr 2 * x) [fc] where x is a signal of the form of (1.31 and g ai is the convolution kernel. Therefore, our 
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results hold immediately if the combined convolution kernel g [fc] = (g ai * g a2 ) [&] meets the definition of 
admissible kernel and the signal’s support K satisfies the associated separation condition. For instance, if g 
and g are both Gaussian kernels with standard deviations of o\ and < 72 , then g is also Gaussian with standard 
deviation of a = \Jo\+ o\ and thus obeys the definition of admissible kernel. 

The main theoretical contribution of this paper is the following theorem, stating that under the separation 
condition, for any k m G K with sufficiently large amplitude c m , there exists a close k m G K. Namely, the 


solution of (2.11 locates a spike near any spike of the underlying signal. In ED, the support detection for 
the super-resolution problem was examined by convex optimization tools. We follow this line of research and 
present the main theorem of this paper, as follows: 


Theorem 2.10. Let y be of the form of (1.1) and let g be an admissible kernel as defined in Definition 2.2 


If the sig nal s upport K satisfies the separation condition of Definition \ 2. /j then the solution x of the convex 
program (2.1) has the following property: For any m such that \c m \ > ^D—5 there exists a k m G K such that 


km k r 1 


< 


8j 2 


5(0 )5 


^6 


The localization result of Theorem |2 .1 0| cannot be derived directly from the robustness result of Theorem 
2.5| For instance, suppose that the amplitude c m and the noise level 6 are of the same order. Then, Theorem 
guarantees that there exists at least one location k m G K obeying 


2.10 


k — k 


< C for some kernel- 


dependent constant C. This localization error cannot be inferred from the robustness result. 


3 Proof of Theorem 2.10 


The main pillar of the proof is the existence of a special superposition of the kernel and its first derivative 
which satisfies some interpolation properties. This function is frequently called the dual certificate. Next, we 
need the following Lemma from |7|: 

Lemma 3.1. Suppose that the set T := {t m } C R. obeys the separation condition 

min It,- — t,| > va, 

for some constant v > 0. Let {u m } G ±1 be a sign sequence and let g(t) be an admissible kernel as defined 
in Definition 2.2 Then, there exists a set of coefficients {a m } and {b m } such that the function 


9(*) = J2 


dm 9 


t-t r 


b m g 


( 1 ) 


t t„ 


satisfying = max t6 f |g(i)| < 1 and interpolates the sequence {u m } on {t m }, namely q ( t m ) = u m for all 

t m G T. Additionally, for all t satisfying \t — t m \ < ea for some t m G T , we have |g(i)| < 1 — ~ 0 *™2 ■ 

We will also need the following result: 

Lemma 3.2. Let the signal support K satisf y the sepa ratio n condition of Definition 2.1 and let K near '■= 
^nearDK, where 7j near and I\ are defined in (2.5) and (2.2), respectively. Then, 

64 ff (0) 7 4 


y: \c n \d? (k n , k^j 

{n: fc„e^„ e ar} 


< 


6 





























where d (k, K) := min^gA' | k m 


k | and 7 and 5 are defined in (2.3) and (1.1), respectively. 


Proof. See Appendix [A] 

and let K m := Z m 
get 


We are now ready to prove Theorem 


2.10 


Fix k m £ K , let Z m be as in (2.4 


Since e < v (the separation constant), Z m C\K = {k m }. Hence by Theorem 2.5 


we 


□ 

n K. 



— 

J2 2%] - x[k\ 

{n: k n ek m } 


fcez m 


< ^2 \x[k] - x[k]\ 

fee z m 


< 


< 


fee z 
16 7 % 


x[k}\ 


By the triangle inequality we then get |c m | — 

empty set, namely, there exists at least one k 
Then, using Lemma |3.2| we directly get 


^jr S < E {„-.k n £K m } |Cn|- If |cm| > ^S, then I<„ 


G K so that 


k n k<n 


< Nsa. Let d = min 


k n eKrr 


is not an 
i — k 

A 'n nj m • 


This concludes the proof. 


d 2 < 


< 


64g(0)7 4 

E{„: k„EK near } \Cn\ 
64g(0)7 4 ^ 

/3 2 (km|-^) ' 


4 Ultrasound Experiments 

The theoretical results and aforementioned algorithms were validated in a series of numerical and in-vitro 
experiments. In the in-vitro experiments, a mechanical phantom with a known geometry was scanned using 
a single-element ultrasound transducer. The same setup was used in order to estimate the PSF for the 
numerical experiments. An unfocused, single-element Panametrix transducer with 0.375” diameter was used 
in all the scans (A3265-SU, Whaltham, MA). This transducer was operated in a transmitter-receiver pulse- 
echo configuration. A single channel arbitrary function generator (8024, Tabor Electronics, Tel-Hanan, Israel) 
was used to drive the transducer with a pulse containing adjustable number of sinusoidal cycles. The function 
generator was operated at 4MHz during all the experiments described below. The received RF-lines were 
sampled at a rate of 100MHz using a 14bit digitizer (CS8427 ,Gage, Lockport, IL). In order to protect the 
digitizer from over-voltage the received RF data was passed through an external diplexer (RDX-6, RITEC, 
Warwick, RI). A diagram of the experimental setup can be found in Figure pi~3| ) .The recorded RF signals 
were loaded into a personal computer and analyzed using dedicated MATLAB@ code (The Mathworks, 
Natic, MA). 
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Figure 4.1: Representative results of numerical simulation experiments: reflectivity estimates are presented 
against true reflectors locations, for two different pulse lengths at SNR = 23dB. (a) Results from a two-cycle 
pulse, (b) results from a three-cycle pulse. 


4.1 Numerical Simulation 


In order to emulate a realistic RF line, the ultrasound pulse used in the numerical simulations was the 
pulse-echo signal reflected from a thick perspex sheet immersed in a water tank. Pulses with two and three 
cycles were transmitted and the reflected echos were received and recorded. It is assumed that the received 


pulse measured this way is a reliable estimation for the convolution kernel in equation (1.2 I. The simulated 
ultrasound signal was produced according to the following steps. First, the simulated reflectivity signal 
was produced. The reflectivity signals were comprised of a series of evenly-spaced reflectors with random 
amplitudes uniformly distributed between ±[5,10]. Different separations enabled the comparison between 
the localization results achieved for different reflectors density. Second, the simulated reflectivity function 
was convolved with the estimated kernel. Different levels of SNR were simulated by adding a Gaussian white 
noise to the resulting signal. 100 simulated signals were produced for each combination of separation and 
SNR conditions. 

Although the RF pulses as modulated signals do not decay monotonously, the demodulated signals and 
their envelops are compatible with the model described above Therefore, the simulated signals were 

demodulated by multiplication with the carrier frequency and application of a low-pass FIR filter. Following 


the demodulation, the optimization problem (2.1) was solved using CVX, a toolbox for specifying and solving 
convex programs M- The value of 5 was selected so that ||77 * g\\ e < S in accordance with equation (12.1| 


and adjusted according to the SNR level. Examples of the original and reconstructed reflectivity functions 
produced using a two-cycle pulse and a three-cycle pulse with SNR of 23 dB are presented in Figures |4Tj a) 
and[4Tjb), respectively. This example demonstrates how the localization of the estimation decreases and the 
locations of the recovered reflectors get further from their original locations as the width of the convolution 
kernel increases. 

In order to quantify the performance of the localization algorithm, the localization error was measured for 

















































































Number of layer 

Layer Thickness (mm) 

Layer Spacing (mm) 

1 

3 (Perspex) 




1.3 

2 

1 (polycarbonate) 




0.9 

3 

1.8 (Perspex) 




1.3 

4 

1 (polycarbonate) 




1.5 

5 

1.8 (Perspex) 




1.3 

6 

3 (Perspex) 



Table 4.1: Dimensions of the Perspex-polycarbonate phantom. The gap between each Perspex/ polycarbonate 
layer was filled with water as the phantom was submerged in a water tank. 


different levels of SNR using two-cycle three-cycle pulses. The localization error was defined as the distance 
between the real location of each reflector and the location of the closest reflector in the solution. The results 
of this analysis can be seen in Figure 4.2 For each separation and SNR level, subplots (a) and (b) present 


the mean localization error of the simulations using two-cycle and three-cycle pulses, respectively. The mean 
localization error of the simulations produced using the two-cycle pulse is relatively small for SNR levels 
greater than 15 dB and increases for small separation intervals. In contrast, the mean localization error of 
the simulations produced using the three-cycle pulse is markedly higher as the reflectors are not sufficiently 
separated. Subplots (c) and (d) present the standard deviation of the localization error for each separation 
and SNR level. When measuring the mean localization error most of the localization error is averaged out, 
indicating that the solutions are clustered around the real locations. The measured standard deviation values 
are almost one order of magnitude higher than these of the mean localization error. Hence, the standard 
deviation of the localization error is a much more suitable criterion for localization loss. 


4.2 In-vitro experiment 

In order to validate the algorithm’s ability to reliably detect the location of reflectors in real ultrasound 
RF scans, a phantom with a known geometry was constructed. Similar to the phantom used in pQ, the 
phantom was constructed from parallel layers of perspex and polycarbonate held together by a rigid frame. 
The dimensions of the phantom in the direction orthogonal to the transmitted beam are much larger than 
each layer’s thickness making it a ID phantom. The thickness of the layers and the distance between them 
were measured in advance and are detailed in Table |4.1| The large differences between the impedance of the 
solid perspex and polycarbonate sheets and the water in which they are submerged produce strong acoustic 
reflectors in each interface. Therefore, only a fraction of the acoustic energy penetrates beyond the first 
few layers. Due to the fast decay of the signal only the first five layers (six reflectors) were relevant to our 
analysis. 

The phantom was design to enable the validation of the localization algorithm under a various of relevant 
conditions. The thickness of the Perspex and polycarbonate layers and the distance between them are of the 
same order as the wavelength of the transmitted pulse. Therefore, the RF-lines can have variable levels of 
echo interference as function of the number of cycles transmitted in each pulse. Additionally, although the 
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Figure 4.2: Quantification of the numerical simulations: In subplots (a) and (b) the mean localization error 
is presented for each separation and SNR level for simulations produced using a two-cycle pulse and a three- 
cycle pulse, respectively. In subplots (c) and (d) the standard deviation of the localization error is presented 
for simulations produced using a two-cycle pulse and a three-cycle pulse, respectively. 
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Figure 4.3: Experimental set-up. (a) a diagram of the experimental setup, (b) the mechanical phantom 
placed in front of the transducer. 


alternating high reflectivity values do not mimic a specific tissue of organ, they enable the validation of the 
localization algorithm under various SNR values. 

The transducer was inserted through one of the walls of the a dedicated water tank filled with deionized 
water at room temperature (see Figure |4~3| (b)). The phantom was placed in the center of the tank and was 
aligned to be perpendicular to the beam by changing its orientation until a maximal echo amplitude was 
achieved. Pulses with one to four cycles were transmitted and the reflected echos were received and recorded. 
The pulses containing different number of transmitted cycles produced signals with slightly different peak 
amplitudes. In order to achieve comparable base-ground SNR level, white Gaussian noise was added to the 
received signal resulting in identical SNR of 41dB. The pulse shape was estimated from the echo reflected 
from the first Perspex layer, which was thick enough to ensure the separation of the first echo from the 
consecutive reflections. Later, different levels of SNR were achieved by adding white Gaussian noise with 
different standard deviation to the resulting signals. 

The RF signals are demodulated as described in Section |4~T| After the demodulation, the optimization 
problem (2.11 was solved using CVX [14]. Complex reflectivity function values were enabled in order to 
facilitate the incorporation of phase delays into the signal model as proposed in [33|. The value of the 
noise level parameter 6 was selected so that ||77*ff|| fl < 5 in accordance with equation (2.1) and adjusted 
according to the SNR level. The reflectivity function estimates were compared against the true location of 
the reflectors used as a reference. The absolute value of the estimated reflectivity function is presented along 
with the location of the reflectors. Since both the attenuation and acoustic reflections gradually reduce the 
energy of the traveling wave, there was no attempt to compare the estimated magnitude of the reflectivity 
function to that of the phantom. 

Two RF-lines received after transmitting a single-cycle pulse and a four-cycle pulse are presented in Figures 


4.4 a) and |4.4[b), respectively. The estimated reflectivity function correlates well with the structure of the 


phantom and the different interfaces originally masked by the interference of the pulses are revealed. When 
wider pulses are transmitted, the kernel-dependent minimal separation condition demands greater separation 
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Figure 4.4: Representative results of phantom experiments: reflectivity estimates are presented against true 
reflectors’ location, for two different pulse lengths at SNR = 23dB. (a) Results from a single cycle pulse (b) 
Results from a four cycle pulse. Each subplot contains the received RF-line, the estimated and the original 
reflectivity functions and the original and reconstructed envelopes of the signal. 


of the strong reflectors. As a result, better localization of the reflectivity estimates are achieved for shorter 
pulses, containing one or two cycles. Since the reflected energy is not compensated after each layer transition, 
the amplitude of the estimated reflectivity function decrease with depth. As the pulse energy decrease with 
depth, the local SNR decreases, resulting in reduced localization of deeper phantom layers. 

In order to quantify the performances of the localization algorithm, the localization error was measured for 
different levels of SNR and for pulses containing different number of cycles. 100 experiments were performed 
for each SNR level and pulse length combination. A representative example of an estimation performed with 
different SNR levels can be found in Figure pi~5| As the SNR decreases the localization error increases. For 
sufficiently large SNR, the support of the solution (of the discrete optimization problem) is tightly clustered 
around the true support (on the continuum). This phenomenon is typical for i\ minimization problems (see 
the analysis on compact domains in nsmsi). Additionally, we present a comparison between our convex 
optimization approach and two widespread algorithms: MUSIC and orthogonal matching pursuit (OMP). 
The convex optimization approach shows competitive results and outperforms these algorithms for low SNR 
values. We stress that we do not aim to make a thorough performance comparison of different algorithms for 
the decomposition of stream of pulses. 
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Figure 4.5: Localization reduction for low SNR levels: reflectivity estimates of the 3rd - 6th reflectors are 
presented against true reflectors location, for single cycle pulse at two different SNR values. Subplot (a) shows 
the results for SNR = 29dB and subplot (b) shows the results for SNR = 17dB. Subplot (c) presents the 
standard deviation of the localization error of (2.11 for different pulses as function of the SNR. Additionally, 
the errors of OMP and MUSIC algorithms are presented for a single cycle pulse. 100 experiments were 
performed for each case. 


It is important to note that the measured SNR reflects the SNR level for the closest reflector. The signal 
resulting from reflectors located further from the transducer has lower SNR. As the noise level increases the 
signal reflected from some of the deeper reflector drops below the noise floor. Our main result in Theorem 
|2.10| do not guarantee the detection and localization of the reflectors in these cases. 

The reduced localization for low SNR levels and long pulses is reflected by higher standard deviation of 
the localization error (see Figure |T~5| c)). When transmitting short pulses the localization of the reflectors 
remains stable even for low SNR levels. The standard deviation of the localization error for each SNR level 
is higher for longer pulses, as predicted. 

Following Corollary [T9] the solution x can include some false reflectors detections of small amplitude, far 
from the real location of the reflectors. In order to quantify the amplitude of these estimation artifacts as a 
function of the SNR level, the mean value of between the first and the second reflector was evaluated 
(see gray region in Figure |4.6[a)). The minimal distance from each reflector used in this analysis was 
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Figure 4.6: False detections: The mean value of non-zero reflectivity estimates located between the first and 
second reflector but far from the true location of these reflectors (gray region in subplot (a)) was evaluated 
in order to quantify the estimation artifacts level (subplot (b)) as a function of the SNR. 


two wavelength. Since the distance between the first and second reflections in the measured (low-resolution) 
signal is reduced for longer pulses, the analysis was performed for the single cycle measurements. A monotonic 
reduction in the mean amplitude of the false detections as a function SNR was measured as presented in 
Figure |4~6|(b). This inverse relation between the SNR and the estimation artifacts corroborates Corollary 2.9 


5 Conclusions 

In this work, we derived theoretical results regarding the localization properties of convex optimization 
technique for decomposition of a stream of pulses. This highlights the appealing properties of the convex 
optimization approach for the decomposition of stream of pulses, as proved and discussed in previous works 

mu- 

From the theoretical perspective, we demonstrated the effectiveness of the separation condition. Fur¬ 
thermore, it is clear that without some separation, no algorithm can decompose a stream of pulses in the 
presence of noise. However, it is not clear yet that a separation condition as presented in Definition |2.1| 
is necessary, and what is the optimal regularity condition for the recovery by tractable algorithms, such as 
convex optimization, algebraic methods and iterative algorithms. 

The experimental results presented in Section [4] corroborate our theoretical results. With sufficient sepa¬ 
ration of scatterers, the localization degree of strong scatterers was shown to be governed by the noise: the 
mean and the standard deviation of the localization error increases with the noise level. In both numerical 
simulations and phantom experiments, for a given SNR level, short pulses produced better localization com¬ 
pared to wider pulses. In addition, the amplitude of artifact scatterers located far from the location of the 
real scatterers is shown to be small and decrease monotonically with the noise level. 
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The results presented in this paper have important implications on studies attempting to perform super¬ 
resolution on signals that can be modeled as streams of pulses. Specifically, the results of ultrasound de- 
convolution algorithms based on signal sparsity/compressibility priors can be analyzed and explained. The 
reliable localization of strong reflectors in ultrasound images has important applications such as the accu¬ 
rate measurement of the intima-media thickness (IMT), the innermost layers of the of an artery, in order 
to facilitate early detection of plaque formation H3 mi- Another possible use could be in non-destructive 
tests. Moreover, the stream of pulses model is a popular model in compressed sensing ultrasound mm 
and compressed sensing radar [5]. Therefore, the presented results have implication on the interpretation of 
the representations of the signals in these applications. 
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paper. 
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Proof of Lemma 


3.2 


Let u m = sign(c m ), where sign(z) := A and q(t) the corresponding interpolating function as defined in 


Lemma 3.1 Denoting q[k] = q ( k/N ), k £ Z we get from Lemma 3.1 for all k £ Z. 
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Then, combining (A.l) and (A.21 into 
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we get 
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where we used the fact that Efcez„ ear 9[&]£[&] = S{m fc m eA'} l c m|- To conclude the proof, we use (2.11 to 
derive 
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This completes the proof. 
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